Chaos in Hamiltonians with a Cosmological Time Dependence 
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This paper summarises a numerical investigation of how the usual manifestations of chaos and reg- 
ularity for flows in time- independent Hamiltonians can be alterred by a systematic time-dependence 
of the form arising naturally in an expanding Universe. If the time-dependence is not too strong, the 
observed behaviour can be understood in an adiabatic approximation. One still infers sharp distinc- 
tions between regular and chaotic behaviour, even though "regular" does not mean "periodic" and 
"chaotic" will not in general imply strictly exponential sensitivity towards small changes in initial 
conditions. However, these distinctions are no longer absolute, as it is possible for a single orbit 
to change from regular to chaotic and/or visa versa. If the time-dependence becomes too strong, 
the distinction between regular and chaotic can disappear so that no sensitive dependence on initial 
conditions is manifest. 



PACS number(s): 98.80.-k, 05.45.Ac, 05.45.Pq 
I. MOTIVATION 

In time-independent Hamiltonian systems sharp, qual- 
itative distinctions can be made between two different 
types of behaviour, namely regular and chaotic. Regu- 
lar orbits are multiply periodic; chaotic orbits are ape- 
riodic. Chaotic orbits exhibit an exponentially sensitive 
dependence on initial conditions, which is manifested by 
the existence of at least one positive Lyapunov exponent; 
regular orbits exhibit at most a power law dependence on 
initial conditions jl]]. This distinction is, moreover, abso- 
lute in the sense that it holds for all times: an orbit that 
starts chaotic will remain chaotic forever; a regular orbit 
remains regular. This distinction is important physically 
because it implies very different sorts of behaviour (al- 
though topological obstructions like cantori S can make 
a chaotic orbit "nearly regular" for very long times). 

An obvious question then is to what extent this dis- 
tinction persists in cosmology where, oftentimes, one is 
confronted with a Hamiltonian manifesting a system- 
atic time-dependence reflecting the expansion of the Uni- 
verse. If the Hamiltonian acquires an explicit secular 
time-dependence, one would expect that periodic orbits 
no longer exist; and one might anticipate further that 
"chaotic" need not imply a strictly exponential depen- 
dence on initial conditions. Given, moreover, that the 
form of the Hamiltonian could change significantly over 
the course of time, one might anticipate the possibility 
that a single orbit could shift in behaviour from "chaotic" 
to "regular" and/or visa versa. In other words, the dis- 
tinction between regular and chaotic need not be abso- 
lute. 

Of especial interest is what happens to the gravita- 
tional iV-body problem for a collection of particles of 
comparable mass in the context of an expanding Uni- 
verse. It is by now well known that, when formu- 
lated for a system of compact support which exhibits no 
systematic expansion, the iV-body problem is chaotic in 
the sense that small changes in initial conditions tend 



to grow exponentially. Largely independent of the de- 
tails (at least for N 3> 1), a small initial perturbation 
typically grows exponentially on a time scale t» ~ R/v 
comparable to the natural crossing time for the system. 
Does this exponential instability persist in the context of 
an expanding Universe, or does the expansion vitiate the 
chaos? 

To better understand various phenomena in the early 
Universe, attention has also focused on the behaviour of 
systems modeled as a small number of interacting low fre- 
quency modes (and, perhaps, an external environment, 
typically visualised as a stochastic bath). The resulting 
description is usually Hamiltonian (albeit possibly per- 
turbed by the external environment), but, because of the 
expansion of the Universe, which induces a systematic 
redshifting of the modes, the Hamiltonian is typically 
time-dependent. Were this time-dependence completely 
ignored, as might be appropriate in flat space, the re- 
sulting solutions would divide naturally into "regular" 
and "chaotic." The obvious question is to what extent 
such appellations continue to make sense when one allows 
properly for an expanding Universe? One might be con- 
cerned that such systems are intrinsically quantum, and 
that there is no such thing as "quantum chaos" . Q How- 
ever, at least in flat space classical distinctions between 
regular and chaotic behaviour are typically manifested in 
the semi-classical behaviour of true quantum systems, so 
that one might anticipate that any diminution of chaos 
associated with an expanding Universe could have im- 
portant implications for such phenomena as decoherence 
and the classical-to-quantum transition. S 

Section II discusses the onset of chaos in time- 
independent Hamiltonian systems as resulting from pa- 
rameteric instability and, by generalising the discussion 
to time-dependent Hamiltonians, makes specific predic- 
tions as to how chaos should be manifested in the context 
of an expanding Universe. Section III confirms and ex- 
tends these predictions with numerical simulations per- 
formed for a simple class of models, namely two- and 
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three-degree-of-freedom generalisations of the dihedral 
potential ||. Section IV concludes with speculations on 
the implications of these results for real systems in the 
early Universe. 



II. CHAOS AND PARAMETRIC INSTABILITY 

To facilitate sensible predictions as to the manifesta- 
tions of chaos in time-dependent Hamiltonian systems, 
it is worth recalling why, in a time- independent Hamil- 
tonian system, chaos implies an exponential dependence 
on initial conditions. 

A time-independent Hamiltonian of the form 



ff(r,p)=p 2 /2 + U(r) 

leads immediately to the evolution equation 
d 2 r a ^ g V ^ 



dt 2 



di 



(i) 



(2) 



Whether or not an orbit generated as a solution to this 
equation is chaotic depends on how the orbit responds 
to an infinitesimal perturbation. In particular, one defi- 
nition of Lyapunov exponents is that they represent the 
average eigenvalues of the linearised stability matrix as 
evaluated along the orbit in an asymptotic t — > oo limit. 
|?J To understand whether or not the orbit is chaotic is 
thus equivalent to understanding the qualitative proper- 
ties of the linearised evolution equation, 



d 2 6r a 



d 2 V 



dr b dr a 



Sr b 



(3) 



r (t) 



At each point along the orbit, this relation can be 
viewed as an eigenvalue equation, with different eigen- 
vectors Sr A satisfying 



d 2 Sr A 
~d£ 2 ~ 



= -n 2 A (t)Sr J 



where Vl 2 A represents some (in general) complicated func- 
tion of time. Formally, one can expand Vl 2 A in a complete 
set of modes, both discrete and continuous, so that 



d 2 6r A 
dt 2 



(5) 



where of course the sum is interpreted as a Stiltjes in- 
tegral. The important question is then what sorts of 
solutions this equation admits. 

Were tt A constant and negative, this corresponding to 
an imaginary frequency, it is clear that the solution to 
eq. (4) would diverge exponentially, so that the unper- 
turbed orbit must be chaotic. This is in fact the essence 
of Hopf 's H proof that (in modern language) a geodesic 
flow on a compact manifold of constant negative curva- 
ture forms a if-flow. However, chaos does not require 



that Q A be negative. Indeed, as stressed, e.g., by Pet- 
tini H in a somewhat more formal setting, eq. (4) can 
admit solutions that grow exponentially even if fl A is 
everywhere positive. 

As a simple example, suppose that, in eq. (5), only 
one time-dependent mode need be considered, so that 
(with an appropriate rescaling) this relation reduces to 
the Matthieu equation flM 



d 2 Sr 



dt 2 



= -(a + (3cos2t)6r J 



(6) 



For positive a and \/3\ < a, fl 2 is always positive but, 
nevertheless, for appropriate values of a and /3 eq. (6) 
exhibits a parametric instability which triggers solutions 
that grow exponentially in time. 

Eq. (6) provides a simple way to understand physi- 
cally why, if one probes a curve of initial conditions in 
the phase space of some Hamiltonian system, one finds 
generically that that curve decomposes into disjoint reg- 
ular and chaotic regions. Moving along the phase space 
curve corresponds to motion through the a — (3 Matthieu 
plane, but a generic curve in this plane will typically 
intersect both stable and unstable regions, i.e., regions 
where Sr A remains bounded and regions where Sr A grows 
exponentially. 

Given this logic, it is easy to see what changes might 
arise if the unperturbed evolution equation is altered to 
allow for the expansion of the Universe. Suppose that 
one is considering a Universe idealised as a spatially flat 
Friedmann cosmology, for which 



ds 2 = -dt 2 + a 2 (t)S ab dx a dx b 



(7) 



The natural starting point for most field theories is then 
an action 



S= / d 4 xa 3 {t) 



(9 t $) - a- 2 5 ab d a <S>d b <S> - U($) 



(4) which leads to mode equations of the form 



dt 2 



a d<t> k k 2 

3 — JT = — 2^ k 
a dt ar 



dV 
d(j>k ' 



(8) 



(9) 



Similarly, when formulated in the average co-moving 
frame, the evolution equation for a particle moving in 
a fixed potential yields a proper peculiar velocity satis- 

fyingf3 



dv a 
~~dt 



a 

—v 

a 



a 



dV(r,a(t)) 
dr a 



(10) 



In either case, eq. (2) has been changed in two signifi- 
cant ways, namely through the introduction of a frame- 
dragging term rx a/a and the explicit time-dependence 
now entering into the right hand side. The frame- 
dragging contribution contribution can always be scaled 
away by a redefinition of the basic field variable, but the 
time-dependence on the right hand side cannot in general 
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be eliminated. It follows that, generically, eq. (2) will be 
replaced by 



d 2 r a 
~dW 



dV(v,a(t)) 
dr a 



(11) 



The simplest case arises when the time-dependence enters 
only as an overall multiplicative factor, so that 



J2 _Q 



df 2 



-R[a(t)} 



dV(r) 
dr a 



(12) 



To the extent that perturbations of a solution to eq. (2) 
can be represented qualitatively by the Matthieu eq. (6), 
it is not unreasonable to suppose that perturbations of 
its time-dependent generalisation (12) can be represented 
by a generalised Matthieu equation of the form 



d 2 5r / 
dt 2 



-R(t)(a + j3cos 2t)Sr J 



(13) 



or, perhaps, some generalisation thereof in which It is 
replaced by some T(t). 

The qualitative character of the solutions to this equa- 
tion are easy to predict theoretically and corroborate nu- 
merically, at least when the time-dependence of R is not 
too strong and this time variability can be treated in an 
adiabatic approximation: There is still a sharp distinc- 
tion between stable and unstable behaviour, but "unsta- 
ble" does not in general correspond to strictly exponen- 
tial growth. For the true Matthieu equation (6), unstable 
solutions do grow exponentially, with 5r A oc exp(wi), but 
allowing for a nontrivial R(t) leads instead to solutions 
of the form 



Sr A 



exp 



(14) 



If, e.g., R(t) oc t p , chaos should correspond to the exis- 
tence of perturbations that grow as 



exp[i 



1+P/21 



(15) 



In other words, the evolution of 5r A will be sub- or super- 
exponential. 

Two other points should also be clear. First and fore- 
most, the distinction between regular and chaotic should 
no longer be absolute. Equation (13) can be interpreted 
as an ordinary Matthieu equation with time-dependent 
"constants" a — R}l 2 ct and (3 = R}l 2 fi. However, the 
fact that these "constants" change in time means that 
the equation satisfied by Sr A is drifting through the a — (3 
plane, so that the solutions can drift into and/or out of 
resonance. In other words, the behaviour of a small per- 
turbation can in principle change from stable to unstable 
and/or visa versa, which corresponds to the possibility of 
transitions between regular and chaotic behaviour. 

The other point is that, if the time-dependence is too 
strong, the adiabatic approximation could fail and the 
expected behaviour could be quite different. In partic- 
ular, for R oc t~ p , with p — » —2, eq. (15) implies that 



the subexponential evolution expected for p somewhat 
smaller than zero will degenerate into a simple power 
law evolution where, seemingly, all hints of chaos are lost. 
An important question would seem to be how negative p 
must become before the distinctions between chaos and 
regularity are completely obliterated. 

Corroboration of the behaviour predicted in eq. (15) 
and a partial answer to this last question are provided 
by the numerical computations described in the following 
section. 



III. NUMERICAL SIMULATIONS 

The experiments described here were performed for 
time-dependent extensions of the potential 

V (x, y, z) = -(x 2 +y 2 + z 2 ) + ^(x 2 + y 2 + z 2 ) 2 
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/ 2 2 , j, 2 2 , 2 2\ 

(y z + bz x + cx y J, 



(16) 



with variable parameters b and c of order unity, which 
constitutes a three-dimensional analogue of the so-called 
dihedral potential 0. For b = c = 1, the case treated 
in greatest detail, this corresponds to a slightly cubed 
Mexican hat potential. At high energies, the potential 
is essentially quartic, although the quadratic couplings 
ensure that there are significant measures of both regular 
and chaotic orbits. At relatively low energies, somewhat 
less than zero, orbits are confined to a three-dimensional 
trough where, qualitatively, they behave like orbits in 
the "stadium problem," |l2| scattering off the walls in a 
fashion that renders them largely chaotic. 

Some of the experiments involved allowing for a time- 
dependence of the form 



V(x,y,z,t) = R(t)V Q (x,y,z) 



(17) 



with R(t) oc t p . Others involved the alternative time- 
dependence 



V(x, y, z, t) = V [R(t)t, R(t)y, R(t)z] 



(18) 



again with R(t) oc t p . 

Individual sets of experiments involved freezing the 
energy at some fixed value, usually E[V ] = 10.0, and 
selecting an ensemble of some 1000 to 5000 initial con- 
ditions. Three-degree-of-freedom initial conditions were 
generated by setting x — 0, z — const, uniformly sam- 
pling the cncrgcticaly allowed regions of the y — z—p y —p z 
hypercube, and then solving for p x — p x (y, z,p y ,p z , E). 
Two-degree-of-freedom initial conditions were generated 
by setting x — z = p z = 0, uniformly sampling the ener- 
getically allowed regions of of y — p v plane, and solving 
for p x = p x {y,Py, E). The computations were started 
at t = 1, t = 10, or t = 100 and ran for a total time 
T = 256 or T = 512, this corresponding in the absence of 
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any time-dependence to ~ 100 — 200 characteristic cross- 
ing times. The orbital data were recorded at intervals 
At = 0.25. An estimate of the largest short time Lya- 
punov exponent |l3f| was computed in the usual way || by 
simultaneously evolving a slightly perturbed orbit (initial 
perturbation Sx = 1.0 x 10~ 8 ) which was periodically 
renormalised at intervals At = 1.0 to keep the perturba- 
tion small. This led for each orbit to a numerical approx- 
imation to the quantity 



y(i) = lim — In 

^ ' «5Z(0)^0 t 



\\sz(t)\\ 

\\SZ(0)\\ 



(19) 



for a phase space perturbation SZ = (\Sr\ 2 + ISp] 2 ). 1 ^ 2 

Since interest focuses on the probability that "chaos" 
does not correspond to a purely exponential growth, the 
data were reinterpreted by partitioning the Lyapunov 
data into bits of length At = 1.0, each of which probed 
the growth of the perturbation during the interval At^ : 

X (A U ) = Xfe + At)fe + At)- X (t^_ (20) 
These bits were recombined to generate partial sums 



At 



At 



In 



\\SZ(U + t ) 
\\SZ(t )\\\ 



(21) 



which capture the net growth of the logarithm of the 
perturbation. These partial sums were then fit to a poly- 
nomial growth law 



£ = A + Bt q 



(22) 



A purely exponential growth would yield q = 1; sub- 
and super-exponential growth correspond, respectively, 
to q < 1 and q > 1. If the perturbation only grows as a 
power law, eq. (22) should not yield a reasonable fit. 

The principal conclusion of the computations is that, 
overall, the adiabatic approximation works very well so 
that, for a broad range of values of p in the potential 
(17), eq. (15) appears satisfied. This is illustrated in Fig. 
1 (a), which exhibits the best fit exponent q as a function 
of p. The data in this plot combine experiments for both 
two- and three-degree-of-freedom chaotic orbits which al- 
lowed for several different values of b and c. Changing 
the values of b and/or c can significantly change the rate 
at which a small initial perturbation grows, thus altering 
the characteristic value of B in eq. (22), but the exponent 
q seems independent of these details. 

As a concrete example, consider the special case of 
two-degree-of- freedom orbits. For p — 0, i.e., no time- 
dependence, initial conditions with E[Vq] ~ 10.0 evolved 
in (17) exhibit a clean split into regular and chaotic, with 
some 72% of the orbits chaotic and the remaining regu- 
lar. Moreover, for this initial energy cantori are relatively 
are unimportant, so that there are few "sticky" chaotic 
orbits trapped near regular islands. N[£], the distribu- 
tion of the final values of £, is strongly bimodal, and it is 



almost always easy to distinguish visually between reg- 
ularity and chaos. For the chaotic orbits one finds (as 
must be the case) excellent agreement with eq. (15). 

As p increases to assume small positive values, the final 
£'s continue to yield a bimodal distibution, N[£), which 
indicates that, in terms of their sensitive dependence on 
initial conditions, the orbits still divide into two distinct 
classes. However, there is a systematic decrease in the 
abundance of regular orbits, so much so that, for p > 0.3, 
there are few if any regular orbits. Moreover, a detailed 
examination of individual orbits show that they can ex- 
hibit abrupt changes in behaviour between chaotic inter- 
vals where £ grows as t 1+ p/ 2 and regular intervals where 
£ exhibits little if any growth. For 0.3 < p < 0.6 (al- 
most) every orbit in the evolved ensembles of initial con- 
ditions seems a member of a single chaotic population 
with q = 1 + p/2. 

However, for larger values of p orbits which began as 
chaotic subsequently exhibit abrupt transitions to regu- 
larity and remain regular for the duration of the integra- 
tion. This behaviour reflects that fact that, at sufficiently 
late times, orbits with p > become trapped in one of 
the four global minima of the potential with Vq(xo, yo) = 
—4/3, located at x = ± -y/4/3, y = i^/4/3, and os- 
cillates in what is essentially an integrable quadratic po- 
tential. Indeed, when the orbit becomes trapped close 
enough to one of these minima, so that V is strictly neg- 
ative, the orbit quickly loses energy until it comes to sit 
almost exactly at rest. This implies that, eventually, £ 
decreases with time, this reflecting the fact that all or- 
bits in the given basin of attraction are driven towards 
the same final point. 

As p decreases from zero to somewhat negative values, 
two more or less distinct populations again appear to 
persist, at least initially, although now the relative abun- 
dance of chaotic orbits decreases rapidly. For relatively 
small values of p, say p < —0.5 or so, £ grows so slowly 
that the chaotic and regular contributions to N\£\ exhibit 
some considerable overlap, and it is no longer easy in ev- 
ery case to distinguish regular from chaotic. For p < —1.0 
the relative measure of chaotic orbits appears to be ex- 
tremely small, and forp < —1.8 it is unclear whether any 
"chaotic" orbits exist at all. To the extent that "chaotic 
orbits" do exist, they are nearly indistinguishable from 
"regular" orbits. 

These general trends are manifested in Figs. 2 (a) - (f), 
which exhibits the final N[£] at t = 266 for integrations 
with, respectively, p = —1.8, p = —0.5, p = 0.0, p = 0.2, 
p = 0.5, and p = 2.5. Further insights into the behaviour 
of small perturbations can be inferred from Figs. 3 (a) - 
(f), which plot £(t) for 150 randomly chosen orbits from 
each of these integrations. 

This behaviour can be contrasted with what obtains 
for two-degree-of-freedom orbits for the same potential 
Vo if one now allows for the time-dependence given by 
eq. (18). Here once again one finds that changing p leads 
to sub- or super-exponential sensitivity, and, as is evident 
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from Fig. 1 (b), that increasing p tends to yield larger 
values of q. Moreover, there is often tangible evidence 
for two distinct types of orbits, seemingly chaotic and 
regular. However, the details change appreciably. 

In this case, one finds that increasing p from zero to 
slightly positive values tends to increase the overall abun- 
dance of regular orbits, and that even those orbits which 
seems chaotic overall often exhibit regular intervals dur- 
ing which £ exhibits little, if any, growth. For values as 
alrge as p = 0.4, it seems that most - albeit not all - of 
the orbits are regular nearly all the time. However, for 
somewhat larger values of p the relative importance of 
chaos again begins to increase, although one continues to 
observe regular intervals. For sufficiently high values of 
p, one finds that, as for the potential (17), all the orbits 
eventually get trapped near one of the four minima of the 
potential, at which point the orbits become completely 
regular and £ becomes negative. 

Alternatively, if one passes from zero to negative val- 
ues of p, one finds that the relative abundance of regular 
orbits rapidly decreases. Indeed, for p < —0.1 or so there 
appear to be virtually no regular orbits. However, the 
chaos is vitiated in the sense that the growth of a small 
perturbation is decidedly subexponential. Indeed, as p 
decreases, one quickly reaches a point where the depen- 
dence on initial conditions is so weak that, even though 
it seems reasonable to term the orbits chaotic, that chaos 
must be viewed as extremely weak. 

This behaviour is exhibited in Figs. 4 (a) - (f), which 
plot N[£] at t = 266 for integrations with, respectively, 
p = -0.6, p = -0.3, p = -0.05, p = 0.5. p = 0.9 and 
p = 1.25. Figs. 3 (a) - (f) plot £(t) for 150 randomly 
chosen orbits from each of these integrations. 

In terms of their sensitive dependence on initial condi- 
tions orbits in these time-independent Hamiltonian sys- 
tems tend to divide relatively cleanly into two distinct 
classes; and, for the case of a time-dependence of the 
form given by eq. (17), those that exhibit such a sensitive 
dependence are well fit overall by the scaling predicted 
by eq. (15). However, to justify designating these orbit 
classes "regular" and "chaotic," it is important to ver- 
ify that these distinctions also coincide with the general 
shapes of the orbits, as manifested visually or through 
the computation of a Fourier transform. 

This is in fact easily done. Orbits that are chaotic 
in terms of their sensitive dependence on initial condi- 
tions tendy systematically to look "more irregular" and 
to have "messier" Fourier spectra than do those which 
do not manifest a sensitive dependence on initial condi- 
tions. For time-independent Hamiltonian systems, one 
can distinguish between regular and chaotic by deter- 
mining the extent to which most of the power is concen- 
trated in a few special frequencies: in a t — > oo limit, the 
power spectrum for a regular orbit will only have sup- 
port at a countable set of frequencies whereas a chaotic 
orbit will typically have power for a continuous range of 
frequencies. However, this is not the case for a time- 
dependent Hamiltonian system! As the energy changes 



in time, the characteristic frequencies of an orbit must 
change so that a long time integration necessarily yields 
broad band power. However, what is true for a regu- 
lar orbit is that this power tends to vary comparatively 
smoothly with frequency. 

This can be understood easily in the context of an adi- 
abatic approximation. At any given instant, it makes 
sense to speak of the two principal frequencies that de- 
fine (say) a two-degree-of-freedom loop orbit; but, as the 
energy changes, the values of these two frequencies will 
change continuously with time. Indeed, the phase space 
can eventually evolve to the point that an orbit that 
starts as a loop becomes unstable and is transformed into 
a chaotic orbit or, perhaps, a different type of regular or- 
bit. 

Representative examples of regular and chaotic orbits 
in time-dependent Hamiltonians are given in Figs. 6 and 
7 which, respectively, exhibit data generated for the po- 
tential (17) with p — —0.4 and the potential (18) with 
p = 0.3. [pTj In each Figure, the left hand panels corre- 
spond to two-degree-of-freedom chaotic orbits; the right 
hand sides correspond to two-dcgrcc-of-frccdom regular 
orbits. The top panels display the orbits in the x — y 
plane. The middle panels exhibit the power spectra, 
|a;(o;)| and |y(w)|. The bottom panels show the total 

IV. PHYSICAL IMPLICATIONS 

At least for the time-dependent potentials described in 
this paper, it seems possible to make sharp distinctions 
between regularity and and chaos, even if these distinc- 
tions are not absolute: as its energy changes, an orbit can 
evolve from regular to chaotic and/or visa versa. How- 
ever, the time-dependence alters the exponential depen- 
dence on initial conditions manifested in the absence of 
any time-dependence, yielding instead a sub- or superex- 
ponential sensitivity. For the models described in Section 
III, p < corresponds physically to an expanding Uni- 
verse; and, as such, the computations corroborate the 
physical expectation that the expansion should suppress 
chaos. Whether or not in the real Universe this expansion 
is strong enough to kill the chaos completely depends on 
the details of the assumed behaviour of both the scale 
factor and the dynamical model. 

But why should one care? Why would the presence 
or absence of chaos matter? Perhaps the most impor- 
tant implication of chaos for the bulk evolution of self- 
gravitating systems is its potential role in violent relax- 
ation, the coarse-grained evolution of a non-dissipative 
system towards a well-mixed state. As formulated origi- 
nally by Lynden-Bell, [^5) violent relaxation is a phase 
mixing process whereby a generic blob of collisionlcss 
matter, be it localised in space or characterised by any 
other phase space distribution, tends to disperse until it 
approaches some equilibrium or steady state. The cru- 
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cial question is one of efficiency. How quickly, and how 
efficiently, will the matter disperse? At a given level of 
resolution, how long must one wait before the matter 
constitutes a reasonable sampling of a near-equilibrium? 

The important point, then, is that the answer to these 
questions depends crucially on whether the flow be reg- 
ular or chaotic. If, in the absence of an expanding 
Universe, the matter executes regular, nonchaotic trajec- 
tories, the approach towards equilibrium will be compar- 
atively inefficient. As probed by a coarse-grained distri- 
bution function and/or lower order phase space moments, 
there is a coarse-grained evolution towards equilibrium, 
but it proceeds only as a power law in time. If, however, 
the matter executes chaotic orbits, this approach will be 
exponential in time, proceeding at a rate that correlates 
with the values of the positive Lyapunov exponents. 

The obvious inference is that to the extent that the 
expansion of the Universe weakens chaos, it should also 
weaken this chaotic mixing. This suggests that phase 
mixing will be a comparatively inefficient mechanism at 
early times, when the dynamics is dominated by the over- 
all expansion, and that it can only begin to play an im- 
portant role at later times, once an overdense region has 
"pinched off" from the overall expansion and begun to 
evolve more or less independently. 
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FIG. 1. (a) Best fit exponent q for a growth rate oc exp(at 9 ) 
for two-degree-of-freedom chaotic orbits evolving in the poten- 
tial V = t p Vo(r) with a = b = 1. The line yields the prediction 
of an adiabatic approximation, (b) The same, except for the 
potential V = V (t p r). 




FIG. 2. (a) Normalised distribution N[£(t = 522)] for 
two-degree-of-freedom orbits with a = b = 1 evolved for the 
interval 10 < t < 266 in the potential V = V (t p r) with 
p = -1.8. (b) The same for p = -0.5. (c) p = 0.0. (d) 
p = 0.2. (e) p = 0.5 (f)p = 2.5. 



FIG. 3. (a) £(t) for 150 representative two-degree-of-free- 
dom orbits evolved in the potential V = t p Vo(r) with 
p = -1.8 (b) The same for p = -0.5. (c) p = 0.0. (d) 
p = 0.2. (e) p= 0.5 (f) p = 2.5. 




FIG. 4. (a) Normalised distribution N[£(t = 522)] for 
two-degree-of-freedom orbits with a = b = 1 evolved for the 
interval 10 < t < 266 in the potential V = V (t p v) with 
p = -0.6. (b) The same for p = -0.3. (c) p = -0.05. (d) 
p = 0.5. (c) p = 0.9 (f) p= 1.25. 



FIG. 5. (a) £(t) for 150 representative two-degree-of-free- 
dom orbits evolved in the potential V = Vo(t p r) with 
p = -0.6 (b) The same for p = -0.3. (c) p = -0.05. (d) 
p = 0.3. (e) p = 0.9 (f) p = 1.25. 




FIG. 6. (a) A chaotic two-degree-of-freedom orbit evolved 
in the potential (17) with p = —0.4. (b) A regular orbit 
evolved with the same value of p. (c) The Fourier transformed 
quantities 10 _3 |x(a;)| and 10 _3 |«/(cj)| (the latter translated 
upwards by 0.5) generated from the orbit in (a), (d) The 
same for the orbit in (b). (e) £(£) for the orbit in (a), (f) The 
same for the orbit in (b). 



FIG. 7. (a) A chaotic two-degree-of-freedom orbit evolved 
in the potential (18) with p — 0.3. (b) A regular orbit evolved 
with the same value of p. (c) The Fourier transformed quan- 
tities and |y(w)| (the latter translated upwards by 40) 
generated from the orbit in (a) . (d) The same for the orbit in 
(b). (e) for the orbit in (a), (f) The same for the orbit 
in (b). 



